Interplay of the Aharonov-Bohm effect and Klein tunneling in graphene 
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We numerically investigate the effect of Klein tunneling on the Aharonov-Bohm oscillations in 
graphene rings using a tight-binding model with nearest-neighbor couplings. In order to introduce 
Klein tunneling into the system, we apply an electrostatic potential to one of the arms of the 
ring, such that this arm together with the two adjacent leads form either a nn'n- or npn-junction 
{n,n: conduction band transport, p: valence band transport). The former case corresponds to 
normal tunneling and the latter case to Klein tunneling. We find that the transmission properties 
strongly depend on the smoothness of the pn-interfaces. In particular, for sharp junctions the 
amplitude profile is symmetric around the charge neutrality point in the gated arm, whereas for 
smooth junctions the Aharonov-Bohm oscillations are strongly suppressed in the Klein tunneling as 
compared to the normal tunneling regime. 

PACS numbers: 72.80.Vp, 73.23.-b, 73.40.Gk, 85.35.Ds 



I. INTRODUCTION 

It is by now common knowledge that graphene has 
peculiar transport properties making the material an in- 
teresting candidate for future applications (for recent re- 
views on graphene see Refs. [iHH)- Ballistic transport in 
graphene has been coined pseudo-difFusive at the Dirac 
poinlFi^~ because it is carried by evanescent modes. This 
gives rise to transport properties that resemble diffusive 
transport in other materials - a prediction that has been 
experimentally confirmed in shot noise measurements 
Phase-coherent transport in disordered graphene is not 
less interesting than ballistic transport in clean graphene. 
The reason is that the honeycomb lattice of graphene 
in combination with different type of scattering mecha- 
nisms yields rather rich localization physicsi^ Depending 
on the magnitude of the so-called intervalley scattering, 
one can either see weak antilocalization or weak localiza- 
tion,^*^ A recent experiment on quantum interference in 
graphene has confirmed this prediction by measuring the 
same sample at different carrier densities and tempera- 
tures which allows to see the transition from localization 
to antilocalizationjii 

Therefore, the combination of ballistic transport at (or 
close to) the Dirac point with quantum interference ef- 
fects suggests itself to contain interesting physics. This is 
our motivation to study how the Aharonov-Bohm effectiS. 
in graphene rings is affected by tuning one of the arms of 
the ring with an external gate through the Dirac point, 
see Fig. [1] for a schematic. We will show below that such 
a setup allows for a clear graphene-specific signature in 
Aharonov-Bohm measurements which seems to be read- 
ily observable. Its physical origin is the quantum inter- 
ference of normal tunneling as well as Klein tunneling 
trajectories through the two arms of the ring. 

Previous theoretical work on graphene ring struc- 
tures includes the investigation of the electronic prop- 
erties of closed single-layer— and bi-layerJd. rings pierced 




FIG. 1: (Color online) Schematic of the graphene ring struc- 
ture (left) and the y-dependence of the on-site gate potential 
V (right) that is applied to the lattice sites on the lower arm 
of the ring. V is non-zero only within the blue box and ex- 
hibits either a smooth (black solid line) or sharp (red dashed 
line) profile along y-direction while being constant along x- 
direction. The width of the arms of the ring is chosen equal 
to the width w of the leads. The shaded area indicates the 
region of non-vanishing homogeneous magnetic field pointing 
out of plane. Different edge disorder configurations are re- 
alized by randomly removing sites within the two regions of 
width 5 between dashed and solid circles. 



by a magnetic fiux as well as transport studies of the 
Aharonov-Bohm effect in clean and disordered graphene 
away from the Dirac point J^ii^ It should be mentioned 
that none of the graphene-specific predictions of Refs.lTsI- 
Tg' have been observed in subsequent Aharonov-Bohm 
oscillation measurements.'^'''-^ Recently, Katsnelson has 
studied the Aharonov-Bohm effect in undoped graphene 
(at the Dirac point) and made graphene-specific predic- 
tions that are complementary to ours and might be ob- 
servable in future experiments.^^ 

The article is organized as follows. In SecHH we intro- 
duce the tight-binding model that we use for the trans- 
port analysis. This section includes a brief description 
of the recursive Green's function formalism with an em- 
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phasis on the peculiarities due to graphene's honeycomb 
lattice. Subsequently, in Sec. Illli we discuss our results 
in the different transport regions which show the inter- 
play of the Aharonov-Bohni effect and Klein tunneling 
in phase-coherent graphene nanostructurcs. Finally, we 
conclude in Sec. IIVI Some technically details are dis- 
cussed in the appendices. 
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II. MODEL 

Our calculation starts with the usual tight-binding 
Hamiltonian for graphene 

n = J2v.\^){^\ + Y.n,\^){J\ (1) 

where the second sum runs over nearest-neighbors and 
Vi = Viji) is an on-site potential that may depend on 
position. The graphene hopping integral tq ~ 3 eV picks 
up a Peierls phase in the presence of a magnetic field 
yielding for the nearest-neighbor coupling element the 
expression 

nj = -To exp A(r) dr^ (2) 

where the line integral is taken along the straight path 
between sites i and j. = h/e is the magnetic flux 
quantum. 

The system under consideration is a ring-shaped struc- 
ture cut out of a graphene sheet, which is attached to two 
crystalline leads also modeled using the graphene lattice 
structure (see Fig. [1]). The structure is subject to a ho- 
mogeneous magnetic field B(r) = (0, 0, B) resulting from 
a vector potential A(r) = {—By,0,0) as well as a gate 
electrode potential Vg located on top of the lower arm 
of the ring. The smoothness of the potential interface 
is controlled via the smoothing width Wg measured from 
the lower edges of the leads: 

V = for y> -w/2, 
V — Vg for y<—w/2 — Ws, 
< V <Vg otherwise, 

taking the origin of coordinates at the center of the ring. 
In our simulations, we used a cosine-shaped smoothing 
profile and chose values Wg = . . . R — 3w/2. 

For a Fermi energy E > 0, together with the adja- 
cent leads this lower arm forms either a nn'n- or npn- 
junction for < iJ and Vg > E, respectively (see Fig. [5] 
for a schematic). Note that the setup exhibits a fiat po- 
tential profile for trajectories along the upper ring arm, 
i.e. a nrm-junction, since there is no gate potential ap- 
plied. This enables a rather large transmission through 
the ring even when the lower ring arm is tuned towards 
the Dirac point, since transport through the upper arm 
always takes places at an energy distance E away from 
the charge neutrality point. 



FIG. 2: Schematic of the influence of the potential profile 
introduced by Vg on the spectrum of the lower arm of the ring. 
The left hand side shows the normal tunneling case (nn'n- 
junction) and the right hand side the Klein tunneling case 
(npn-junction). The location of the van Hove singularities at 
±To is also schematically depicted. 

We derive the transmission function through the 
ring from the scattering 5-matrix using the Landauer- 
Biittiker formalism for elastic transport at zero temper- 
ature assuming complete phase-coherence. The semi- 
infinite left and right leads are described through their 
respective real-space surface Green's functions (SGF) 
and qr. The coupling of the leads via the Hamilto- 
nian H describing the ring structure is expressed through 
Dyson's equation. In order to obtain the effective cou- 
pling between the leads, we apply a recursive Green's 
function (RGF) technique. The 5'-matrix is then ob- 
tained from the system's Green's function using the 
Fisher-Lee relationj2£ 



A. Calculation of the lead's surface Green's 
function 

For the calculation of the lead's surface Green's func- 
tions, we essentially follow Ref. [2l]. An isolated infi- 
nite lead is described as quasi-one-dimensional periodic 
arrangement of identical unit cells. Each unit cell is 
described by an intra-cell matrix Hq whose dimension 
equals the number M of atomic sites within the unit cell. 
Neighboring unit cells are coupled via an inter-cell ma- 
trix Hi , such that for each cell a Schrodinger equation of 
the form 

Hli,,_i+Ho^, + Hi^,+i^Ei,,, zeZ (3) 

holds, where ipz is a M-dimensional vector describing the 
atomic sites within a particular unit cell. Since we con- 
sider the zero temperature regime, transport takes place 
at the Fermi energy E, measured relative to the charge 
neutrality point in the leads. The approach outlined in 
Ref. [21I then requires the inversion of the inter-cell ma- 
trix Hi. However this matrix is singular in the case of 
a graphene lattice. To overcome this problem, we will 
use an effective description for the leads as outlined in 
App. El This description is based on the fact, that in 
many cases not all of the atomic sites contained in one 
cell couple to the adjacent cell. One can then conve- 
niently divide the unit cells into even- and odd-indexed 



3 



subcells and eliminate the latter from the calculation by 
writing down the effective coupling between even-indexed 
subcells. This yields an expression of the form 



(4) 



with invertible coupling matrices Hi. As will be shown 
later, this effective description also has the additional 
benefit of increased performance of the recursive Green's 
function scheme. Note that since the approach does not 
depend on details of the matrices Hq and Hi it may also 
be applied to other lattice structures that exhibit non- 
invertible coupling matrices or that enable a potential 
performance gain. 

The isolated leads' SGFs and gi^ arc thus obtained 
according to expressions (2.16) and (2.17) of Ref. 2l!, re- 
spectively, using the effective description ^ for the leads. 



B. Connecting leads and conductor 

In order to obtain the S'-matrix, we need to compute 
the Green's function 
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describing the effective coupling between the surface sites 
of the two leads. To this end, we apply a variant of 
the commonly-used RGF schemei^"— As described in 
App. [B] in an effective description, the leads have to 
couple to the conductor through an additional contact 
slice, which has to be taken into account in contrast to 
conventional RGF algorithms. 

Similar to the leads, the sample itself is also divided 
into slices. The width Ax of the slices is chosen as small 
as possible, since the algorithm scales only linear with the 
length of the sample (along x-direction) but up to third 
power with the width (along y-direction). We therefore 
choose /S.X = ax/2 for zigzag leads and Aa; — ax/^ for 
armchair leads, where ax is the lattice constant of the 
honeycomb lattice along the a:-direction, which is ax — 
coa/S in the zigzag case and ax = 3ao in the armchair 
case, ao = 0.142 nm being the nearest-neighbor distance 
in graphene. 

The applied RGF procedure then has the following 
structure: 

1. Connect the contact slices to the initially isolated 
leads. 

2. Connect to slice 1 of the conductor. 

3. Set n=l. 

4. Connect to slice n + 1 of the conductor and elimi- 
nate slice n from the description. 

5. Increase n by one and repeat the previous step until 
all slices are connected. 



In each step it is sufficient to update the contact slice 
Green's functions. The Green's function ([5]) of the fully 
coupled leads is updated once when all slices are con- 
nected42, The scheme of a particular step in the recur- 
sion is as follows: Interpreting the coupling to slice n as 
perturbation to the system of coupled slices up to slice 
n — 1, the Green's function describing the coupled system 



9c {n) gcx{n) gcoin) 
9{n) = I 9xc{n) gx{n) gxD{n) 
9Dc{n) goxin) gain) 

is obtained from the unperturbed Green's function 

gc{n-l) gcoin -1) 
go{n)=\ gncin-l) goin-l) 

(i;-ij("))-i 

and the coupling matrix (the perturbation) 

/ To,„ 

U = [ T„_i,„ 

V(To,„)^ (r„_i,„)^ 

via Dyson's equation 

9{n) = 9o{n) + go{n)U g{n). 



(6) 



(7) 



(8) 



(9) 



Here, the index C refers to the contact slices, whereas the 
index D refers to the last connected slice of the sample, 
i?^") is the layer-local Hamiltonian of slice n of the sam- 
ple, and the matrices Tm,n describe the coupling of the 
contact slices and slice n — 1 to slice n of the sample. For 
a simple two-terminal setup as we consider here, To_„ is 
non-zero only for n € {1, A^}, A^ being the total number 
of slices of the sample. 

Note that the matrix elements containing an index X 
in do not have to be calculated since they do not 
appear in the next step of the recursion. Further note 
that we do not need to add an inconvenient infinitesimal 
imaginary part to the energy in expression ([7]) since it 
would anyway be absorbed by finite imaginary terms in- 
troduced by the fact that we deal with an open quantum 
system. 



C. Calculation of the scattering matrix 

The linear conductance of the system is obtained using 
the Landauer formula 



g^g^ Tr(tTi) 



(10) 



where g^ = 2e^//i. The factor 2 accounts for spin de- 
generacy and t is the transmission matrix element of the 
scattering matrix 
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which itself is a matrix whose elements are the trans- 
mission amplitudes for scattering between the different 
transverse modes in the two leads. The ^-matrix can be 
written in terms of the Green's function Gij in Eq. ([5]) 
by means of the Fisher-Lee relation^° 

= ^\ ■ 1 + GyV) c^i ■ ^J^i. (12) 

In the latter equation, the transverse eigenvectors <j)i of 
in-moving states in lead j and the duals 4>\^ of the out- 
moving states in lead i as well as the corresponding group 
velocities Vh^i and the matrix V are defined and calculated 
as described in Ref. 111. 



III. RESULTS 

In the following we present transmission properties for 
a ring with R/ao = 300 and w/ao = 60. Edge disorder 
is applied to the ring by randomly removing sites within 
a width 5 from the inner and outer edges of the ring, re- 
spectively (see Fig. [I]). We choose 5/aQ = 1.5 in order to 
keep the edge of the ring as smooth as possible while still 
allowing for different edge disorder configurations. As 
depicted in Fig. [21 Fermi energy i? e ]0, . . . , to[ and gate 
potential Vg € . . . ,2 E] are chosen such that transport 
always takes place in between the van Hove singularities 
located at i? = ±to where the density of states diverges 
in the tight binding model of graphene. 

In Fig. [3l we plot the magnetoconductance at Fermi 
energy E/tq = 0.5 and zero gate voltage [Vg = 0) 
for a particular ring realization, showing pronounced 
Aharonov-Bohm oscillations on top of a low frequency 
background. The background signal results from uni- 
versal conductance fluctuations (UCF) which is typical 
for phase-coherent mesoscopic devices. The behavior is 
in agreement with the observations made in Ref. \\K 
where the authors investigate an even wider magnetic 
field range up to the quantum Hall regime. In Fig. |4l 
we also show the corresponding frequency spectrum ob- 
tained from a Fourier transform of the magnetoconduc- 
tance. The contributions to the Aharonov-Bohm oscil- 
lations are centered around {/S.B e / h)~^ ~ 2.3 • 10^. 
Using TT • AS = h/e, this frequency corresponds to a 
mean radius R/a^ ^ 270 of interfering electron trajec- 
tories, which perfectly lies within the boundaries of the 
ring. 

In Fig. [5l we show the same plot for E/tq = 0.1. The 
oscillations diminish at Ba^e/h ^ 6 • 10^^. This field 
strength marks the onset of the quantum Hall regime, 
where the cyclotron diameter becomes comparable to the 
width of the arms of the ring; an estimate of the graphene 
cyclotron diameter dc = 2E /vpeB^ taking the Fermi ve- 
locity at the Dirac point in graphene, vp = 37rroao//i, 
yields dd ^ 40, a value of same order of magnitude as 
the width w/ao — 60. 

By applying a gate voltage Vg > to one of the ring 
arms, the magnitude of the Aharonov-Bohm oscillations 




FIG. 3: (Color online) Magnetoconductance of a ring with 
R/ao = 300, w/ao ~ 60 at energy E/tq — 0.5 and zero gate 
voltage, showing clear Aharonov-Bohm oscillations on top of 
a background due to universal conductance fluctuations. 
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FIG. 4: (Color online) Frequency spectrum corresponding to 
Fig. [3l obtained from the Fourier transform Q of the magne- 
toconductance Q. Besides the low frequency background and 
the fundamental oscillation component, the second harmonic 
is also slightly visible in the spectrum. The dashed line in- 
dicates the frequency limit of the high pass frequency filter 
used for background subtraction. 



may be modified. A convenient measure of the oscilla- 
tion magnitude is the root mean square (RMS) ampli- 
tude of the signal. Prior to the RMS analysis, the UCF 
background has to be removed from the signal. This 
is achieved by applying a high pass frequency filter to 
the Fourier transform of the magnetoconductance data, 
as indicated in Fig. The retained, unbiased signal is 
squared, and the root of the average over the squared 
signal is defined as the RMS amplitude AQums- 

In Fig. ini we show the dependence of the RMS os- 
cillation amplitude AQfjMS on the gate voltage Vg for 
different smoothing widths Ws (see Fig. [T]) at energy 
E/tq — 0.5, where the average is taken over the full 
range B — . . . 10~^ $o/ao- Increasing the gate volt- 
age from zero towards the neutrality point Vg = E not 
only leads to increased potential scattering but also to a 



5 




2 4 6 8 10 



FIG. 5: (Color online) Magnetoconductance of a ring with 
R/ao = 300, w/ao — 60 at energy E/to — 0.1 and zero gate 
voltage, showing the onset of the quantum Ifall regime. Note: 
The conductance is still finite near zero magnetic field, which 
is not visible on this scale. 
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FIG. 6: (Color online) RMS analysis for the setup 
used in Fig. [3] for different smoothing widths Ws/ao € 
{0, 21, 52.5, 105, 210}. Each data point results from an aver- 
age over five realizations of edge disorder. The corresponding 
standard deviations lie between 0.005 Go and 0.015 Go but are 
suppressed for better visibility. For better clarity, the spec- 
trum schematics (see Fig. [2]) are also included. 
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FIG. 7: (Color online) Magnetoconductance of a ring with 
R/ao = 300, w/ao — 60 tuned close to the neutrality point 
in the lower arm with a smooth potential interface (ws/ao = 
210). The plots (a)-(d) constitute a representative selection 
and are obtained by variation of the microscopic edge disorder 
configuration (which is always chosen randomly) for values 
Vg/E^ 1.00 ±0.01. 



pends on the smoothness of the gate potential. In case 
of a smooth potential, the partial waves in the lower arm 
have to tunnel through a finite region of low density of 
states, where V ^ E (see Fig. [T|), in order to interfere 
with the partial waves traversing the upper arm. The 
lower arm becomes increasingly penetrable as this region 
gets narrower, until it gets transparent in case of a sharp 
potential. This reflects the usual behavior of Klein tun- 
neling phenomena, where the probability for tunneling 
through a pn-junction depends on the smoothness of the 
pn-interfacei2^ 

The described behavior of the RMS amplitude is ro- 
bust over the whole energy range under consideration, 
except for an increasing uncertainty at lower values for 
the Fermi energy. Although all results are presented for 
zigzag boundary conditions in the leads, the effects are 
independent of a change of orientation of the graphene 
lattice to armchair boundaries in the leads. 



reduction in the number of accessible propagating states 
in the lower arm of the ring. As can be seen in Fig. [51 the 
oscillation amplitude diminishes and reaches a minimum 
value at the neutrality point. Note that, since the trans- 
mission through the upper ring arm is not at all affected 
by a gate potential, the overall conductance itself is only 
slightly changed to fluctuate around 2.5^0 (see Fig. [7]) 
as compared to values around 3.4 t/o in the case of zero 
gate potential on the lower ring arm (see Fig. [S]). 

For Vg < E, the decay of the RMS amplitude towards 
the neutrality point does not depend on the details of the 
gate potential interface. However, in the regime of Klein 
tunneling, Vg > E, the oscillation behavior strongly de- 



Before we conclude, we mention here an additional ob- 
servation concerning the dependence of the magnitude 
of the Aharonov-Bohm oscillations on the magnetic field 
strength B, when the lower arm of the ring is tuned 
near the neutrality point (see Fig. [7]). It seems that in 
this regime the oscillation magnitude is in general signif- 
icantly lower for low field strength, compared to the os- 
cillations at higher field strength. This is indeed the case 
for most of the ring realizations we investigated, though 
not for all of them (see Fig. El^d)). The reason for such a 
behavior is so far not understood. Since the increase in 
oscillation magnitude cannot be related to any particu- 
lar length scale a connection to the quantum Hall effect 
seems unlikely. 
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IV. CONCLUSIONS 

In summary, we have numerically analyzed transport 
through graphene ring structures in the presence of a per- 
pendicular magnetic field based on the recursive Green's 
function formalism. In order to understand the physics 
of the interplay of the Aharonov-Bohm effect and Klein 
tunneling in graphene, we have looked at the influence of 
a local gate over one of the arms of the ring on magne- 
totransport. By varying the gate voltage, we have been 
able to tune this arm from the n-type to the p-type trans- 
port regime via the Dirac point. The analysis of the root 
mean square amplitude of the Aharonov-Bohm oscilla- 
tions clearly shows that the p-type signal is smaller for 
smooth pn-junctions in the ring arm and can recover the 
full n-type value only for very sharp pn-junctions. Our 
predictions nicely complement the analysis of Ref. [l^ 
where both arms of the ring are assumed to be tuned 
to the charge neutrality point. This might lead to the 
first observation of the Aharonov-Bohm effect caused by 
transport through evanescent modes or a combination of 
propagating modes in one arm and evanescent modes in 
the other arm of the ring. 
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Appendix A: Effective description of infinite leads 

The leads may be described in an effective manner if 
some of the atomic sites contained in cell z do not couple 
to the adjacent cell z -I- 1. One can then conveniently 
divide the unit cell into two subcells according to 

resulting in a double-periodic structure as shown in 
Fig. [51 where ft-o only describes atomic sites that directly 
couple to the next cell. Writing i/iz — then 
yields the effective description of the lead ^ where 

Hi = toi {E — hi)^^ tio 
Ho = Oio + ho + 6oi 

eol = toi (E - hi)-Hl^ (A2) 

Note that simple single-periodic structures are included 
in the description via Hq = hi, Iqi — tiQ. 

In the case of a graphene lead with zigzag edges in zero 
magnetic field, the non- vanishing matrix elements of the 



hi ho hi ho hi ho 




z — 1 z z + 1 



FIG. 8: Division of a zigzag (left) and armchair (right) lead's 
unit cell. In the zigzag case, the first of the two slices in which 
the unit cell is divided is described by hi , and the second slice 
that contains those sites which couple to the next unit cell to 
the right is described by ho- Analogously, in the armchair 
case, the first three slices of a unit cell are described by hi, 
whereas only the fourth slice is described by Hq. 

(M/2 X M/2)-matrices toi,tio and hk, k G {0, 1} may be 
written in a compact form, 

(^fe)ii — Vo 

{hk)2i+k-l,2i+k = —To = (hk)2i+k,2i+k-l 

{toi)ii = -To = {tio)u (A3) 

where we assumed tq S M and a constant on-site matrix 
element Vo G M. Of course, the roles of ho and hi are 
interchangeable, dependent on the details of the lead's 
surface. 

Obviously, the effective coupling matrix Hi = Tq (E — 
hi)~^ is invertible, in contrast to the original (M x M) 
coupling matrix Hi in the non-effective description, 
which is clearly non-invertible. Additionally, the size M 
of the matrices describing the effective lead is smaller by a 
factor of two compared to the original description, yield- 
ing increased performance, since the computational effort 
of the applied RGF scheme roughly scales as 0{M^). The 
performance gain is even larger in the case of graphene 
leads with armchair edges, since then only 1/4 of the 
atomic sites within a particular unit cell couple to the 
next cell. The derivation is similar to the zigzag case 
and will not be repeated here. 

One point to note is that Hi is not well defined 
for energy values matching the eigenvalues of hi, E = 
Vo, Vo ±To. E = Vo corresponds to the charge neutrahty 
point in graphene, where the density of states vanishes 
and hence there are no propagating modes in the lead 
that may contribute to the current. The singularity at 
E — Vo i To corresponds to the van- Hove singularity in 
graphene where the density of states diverges. 

For a more formal treatment of this regularization pro- 
cedure see Refs. [2711281 . The subdivision of the graphene 
lead's unit cell as depicted in Fig. [5] has also been previ- 
ously used in Ref. [23. 
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FIG. 9: Comparison of actual (a) and effective (b) leads, to- 
gether with their coupling to a scatterer Hs- Effective leads 
couple to the scatterer through an additional contact slice, 
described by He- 



Appendix B: The lead— sample interface 

For a semi-infinite lead that terminates at z = (see 
Fig. in]) such that the wave function on the surface of the 
lead is described by V'o in the effective description (|4]), 
we have to account for the boundary conditions. For a 
left lead, Eq. (U) still holds for z < and we additionally 



have the equations 



ho -^0 + Til's 



(Bl) 
(B2) 



where V'S and Hs describe the part of the sample that 
couples to the lead, and T accounts for this coupling. 
By eliminating as we have done before, (|Bip may be 
replaced by 



EiPo = Hi V'-2 + Hcipo+Tips 



(B3) 



where He = ho + 9io describes a contact slice, through 
which the coupling of the lead to the sample takes place 
in the effective description, as shown schematically in 
Fig.m 
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